AD-755  401 


INVESTIGATION  OF 
PHYSICS  BY  RADAR 


A 


CLOUD  AND  PRECIPITATION 


J.  S.  Marshall 


McGill  University 


_ J 

- ^ 

Prepared  for: 

Air  Force  Cambridge  Research  Laboratories 
1  March  1972 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Roadr  Springfield  Va.  22151 


AD75540I 


'***--  -*  ^<T,  ■¥>- 


r^amrr^tmrci  ■s.-sty  err* ;  -r?*°  tCT  *-r  ~ 


AFCRL-TR-73-0020 

INVESTIGATION  OF  CLOUD  AND  PRECIPITATION 
PHYSICS  BY  RADAR 
by 

J.  S.  Marshall 

Stormy  Weather  Group 
McGii.1  University 
Montreal,  Canada 


Contract  No.  FI9628-69-C-OIO7 
Project  No.  8620 
Task  No.  862004 

Work  Unit  No.  86200401 


FINAL  REPORT 

1  December  1965  -  28  February  1  72 


March  1972 


Contract  Monitor:  Albert  C.  Chmela 
Meteorology  Laboratory 


c 


Approved  for  public  release;  distribution  unlimited 


Prepared 

for 

AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
BEDFORD,  MASSACHUSETTS  01730 

JUproduc*jd  by 

NATIONAL  TECHNICAL 
INFORAAATION  SERVICE 

U  S  D«po*tm«nt  of  Comm*rc# 

Spnngfi^id  VA  22151 


I 


Best 

Available 

Copy 


1 

i 


Qualified  requestors  may  obtain  additional  copies  from 
the  Defense  Documentation  Center.  All  others  should 
.apply  to  the  National  Technical  Information  Service. 


/ 

t 

! 


* 


v*wmwapaawaCS 


3'2CF*"?*r3?H* 


Unclassified  '*«^**^ 

Security  Classification 

"  ""~  ,  * '  DOCUMENT  CONTROL  DATA  ■  R&D  " 

(Security  classification  of  title,  body  of  abstract  and  indexing' annotation  must  be  entered  when  the  overall  report  is  classified ) 
t.  ORIGINATING  ACTIVITY  (Corporate  author)  ~  ’  in.  REPORT  SECURITY  CLASSIFIcVnoir 

McGill  University  Unclassified 

Stormy  Weather  Group  zb.  group 

Montreal,  Quebec,  Canada 

________  — —  ~ ~  .  .  . 

INVESTIGATION  OF  CLOUD  AND  PRECIPITATION  PHYSICS  BY  RADAR 

4.  OCSCR1PTIVC  NOTES  (Type  of  report  and  inclusive  dates) 

Scientific. _ Final.  1  December  1965  -  28  February  1972  _ 

5.  AUTHOHtS)  (First  none,  middle  initial,  last  name) 

J.  S.  Marshall 


6.  REPORT  DATE 

March  1972 _ 

0O.  CONTRACT  OR  GRANT  NO. 

F19628- 69- C-010T 

&  PROJECT.  TASK.  WORK  UNIT  NOS. 

8620-04-01 

C.  000  ELEMENT  6ll02F 

d.  OOO  5U0ELEMENT  681310 


70.  TOTAL  NO.  OF  PACES  7 tk  NO.  OF  REFS 

34-  I  18 

90.  ORIGINATOR'S  REPORT  NUMBER/SJ 


9b.  other  REPORT  nc(SJ  (Any  other  numbers  that  may  be 
assigned  this  report) 


AFCRL-TR-73-0020 


10.  oistriuution  statement 


A  -  Approved  for  public  release ;  distribution  unlimited. 


ti.  supplementary  notes 


tZ  SPONSORING  MILITARY  ACTIVITY 


TECH,  OTHER 


Air  Force  Cambridge  Research 
Laboratories  (LY) 

L.  G.  Hanscrm  Field 
Bedford,  Massachusetts  01730 


II.  ABSTRACT 


'  Significant  developments  in  instrumentation  directed  toward  the  optimum 
use  of  radar  in  investigations  of  cloud  and  precipitation  physics  './ere 
accomplished.  A  new  technique  in  film  display  of  radar  information  provides 
intensity  maps  in  conformal  map  form.  This  display,  called  AZLOR,  compli¬ 
ments  the  standard  PPI  display  wiffi  detail  and  conformity  to  the  true  pattern 
in  small  areas.  The  technique  of  §>ulse  compression  was  devised  and  temployed  in 
selected  weather  situations.  This  technique  has  the  advantage  of  permitting 
increase  in  sensitivity  without  compromising  the  resolution  of  the  radar; > A 
body  of  microwave  attenuation  statistics  was  accumulated  from  AZLOR  records  of 
severe  storms.  Data  were  accumulated  and  analyzed  to  relate  precipitation 
seen  by  radar  to  cloud  seen  by  photography.  Radio- located  lightning  flashes 
were  related  to  precipitation  by  radar.  A  new  technique  was  explored  to 
extrapolate  radar  precipitation  measurements  to  the  surface  for  comparison 
with  rain  gauges.  Other  research  studies  include;  short  term  forecasting 
using  Height  Plus  Plan  Indicator  records;  peak  reading  as  signal  processing 
technique;  interrelationship  of  showers  an<2  contiguous  rain;  and  measurement 
of  snowfall  radar « 
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The  Stormy  Weather  Research  Group 


The  group's  activities  have  always  been  centred  on  the  McGill  Weather 
Radar.  Since  1956,  however,  the  group's  attention  has  been  shared,  and  is  now 
shared  about  equally,  between  the  McGill  Radar  and  the  Alberta  Hail  Studies. 
Members  of  the  group,  and  their  principal  involvement,  are  as  follows : 


McGill  Radar 
J.S.  Marshall*" 
E.J.  Stansbury 
E.  Ballantyne 
G.L.  Austin 
G.  Drufuca 


Alberta  Hail 

W.  Hits chf eld 

(J.S.  Marshall)-*" 
•f 

R.R.  Rogers 
H.  Leighton 
N.  Cherry 


While  the  two  columns  indicate  a  clear  division  of  responsibilit?  os  (with 
"management"  indicated  by  asterisks),  all  members  of  the  group  and  their 
graduate  students  have  strong  bonds  of  common  interest. 

Prof.  Marshall  now  has  very  nearly  full  time  available  for  research, 
primarily  radar  display  development  and  secondarily  storm  dynamics. 


Dean  Stansbury  concerns  himself  with  radio  location  of  lightning, 
with  a  view  to  relating  lightning  activity  to  the  radar  record  of  precipi¬ 
tation.  As  Dean  of  Science,  he  has  less  than  the  average  professor's  time 
to  devote  to  research. 


Mr.  Ballantyne  is  full-time  engineer-in-charge  at  the  McGill  Radar 
Weather  Observatory,  where  he  has  a  staff  of  one  engineering  assistant 
(Mr.  Abnash  Singh)  and  two  senior  radar  technicians  (Messrs.  Martin  Claaesen 
and  Philipp  Seidenfuss).  He  shared  the  task  of  establishing  the  radar  with 
Dr.  R.W.  Fetter,  who  is  now  with  David  Atlas  at  Chicago.  Display  develop¬ 
ment  is  shared  between  Marshall  and  Ballantyne j  the  latter  did  all  the  work, 
in  that  he  developed  the  circuits  and  with  his  staff  built  the  equipment . 

Dr.  G.L.  Austin  is  assistant  professor  of  physics  and  has  about  half 
his  time  available  for  research.  A  native  of  England  and  a  Cambridge  graduate, 
he  had  experience  in  ionospheric  research  in  New  Zealand  before  coming  to 
McGill  and  redirecting  his  interest  to  the  troposphere.  Apart  from  his 
involvement  in  the  radar,  he  has  taken  over  our  cloud  photogrammetry. 
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Dr.  G.  Drufuca  is  an  electrical  engineer  from  Milan  who  came  to  us 
as  a  research  associate  in  1970.  His  primary  involvement  is  in  obtaining 
from  radar  records  at  a  non-attenuating  wavelength  statistics  on  the 
attenuation  that  would  be  experienced  on  a  higher-frequency  microwave  communi¬ 
cations  link. 

Turning  to  group  members  whose  predominant  involvement  is  in  hail 
research:  Dean  Hitschfeld  is  a  world  leader  in  hall  research.  As  Dean  of 
Graduate  Studies  and  Research,  he  has  less  than  one  quarter  of  his  time 
available  for  research;  he  has  turned  over  the  management  of  our  Alberta  Hail 
undertakings  to  Rogers  and  Marshall,  so  that  the  available  time  will  go  to 
genuine  research.  Prof.  Rogers  was,  like  Dr.  Drufuca,  interested  in  attenua¬ 
tion  statistics,  but  has  been  attracted  to  hail  research  by  the  availability 
in  Alberta  of  an  antenna  uniquely  equipped  for  polarization  experiments. 

Dr.  Leighton  is  a  teaching  fellow  in  the  process  of  converting  himself  from 
nuclear  physicist  to  meteorologist.  Dr.  Cherry,  another  teaching  fellow,  is 
a  young  New  Zealand  physicist  who  studied  lee  waves  in  New  Zealand  and  who 
is  well  suited  to  move  in  rapidly  on  the  problems  of  hail  research. 

The  group  employs  full  time  a  photographic  technician,  a  draughtsman 
and  two  secretaries. 


The  McGill  Radar  Weather  Observatory 

Capital  investment:  a  grant  from  the  National  Research  Council  of 
Canada  paid  for  the  30-foot  antenna  and  pedestal,  radome  and  tower.  It  also 
paid  for  the  modification  of  an  8.6-cm  magnetron  radar,  but  this  was  dis¬ 
allowed  by  the  licensing  authority.  In  this  situation,  efforts  to  obtain  a 
radar  with  transmitting  klystrons  were  renewed,  and  the  F PS-18  was  obtained 
on  loan  through  AFCRL.  The  observatory  building  was  paid  for  by  a  private 
and  anonymous  benefaction  to  McGill  University. 

Sources  of  continuing  support,  additional  to  AFCRL  in  loaning  us  the 
FPS-18,  are  the  National  Research  Council  of  Canada,  the  Atmospheric  Environ¬ 
ment  Service  of  Environment  Canada,  and  the  Communications  Research  Centre  of 
the  Canadian  Department  of  Communications. 


I  .  4 


—  SCIENTIFIC  ATLANTA 
1000  ft.  ANTENNA  RANGE 

O  AT  STE  ANNE  IN  RAOOME 
3  MILE  RANGE  TO  ILE  PERROT 


ELEVATION 


McGill  Weather  Radar  Observatory 
antenna  beam  pattern,  one  way. 
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HEIGHT 


Antenna,  Radome,  Beam  Pattern 


The  antenna  pedestal  is  mounted  on  a  platform  85  ft  above  grade; 
antenna  and  pedestal  are  protected  by  a  geodesic  radome.  With  a  reflector 
of  diameter  30  ft,  the  antenna  gives  a  beam  of  angular  extent  0.8°  in  azimuth 
by  1.0°  in  elevation.  These  dimensions  are  just  to  level  6  dB  (radar)  below 
peak;  for  weather  radar  it  should  be  borne  in  mind  that  width  to  -25  dB  is 
twice  as  great.  Sidelobes  are  below  -40  dB  and  generally  better  than 
-50  dB  (radar). 
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FASE  (Fast  azimuth  slow  elevation)  Antenna  Program 


Since  the  introduction  of  CA^PI  twenty  years  ago,  the  antenna  has 

been  programmed  to  rotate  steadily  in  azimuth,  with  successive  rotations  at 

progressively  higher  angles  of  elevation.  The  present  rate  in  azimuth  is 

six  revolutions  per  minute,  or  tt/5  rad  s-1.  The  low  figure  of  60  s-1  is  used 

for  the  pulse  repetition  frequency,  so  that  the  axes  of  successive  beams  are 

separated  by  rr/300  radian,  or  half  a  beamwidth  at  level  15  dB  below  peak. 

The  elevation  angle  is  held  constant  through  one  rotation  in  azimuth, 

then  raised  incrementally  during  a  quiet  sector  of  2rr/l5  radians  of  azimuth, 

» 

then  again  held  constant  for  one  rotation. 

We  are  now  using  a  GP  (geometric  progression)  program:  the  values  of 
tem  0  (where  0  is  elevation)  increase  in  geometric  progression,  up  to  a  maxi¬ 
mum  of  about  35°  for  0.  The  whole  program  takes  5  minutes,  and  is  repeated 
every  5  minutes. 

The  base  line  from  which  0  is  measured  is  depressed  by  0.5°  below 

the  horizon.  This  depressed  base  line,  when  rotated  about  a  vertical  axis, 

.  * 

sweeps  out  a  conical  approximation  to  the  curved  surface  of  the  earth. 
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POWER  LEVEL 


B 

0 


Power  returned  from  an  extended  target  for  a  pulse  of  constant 
duration  compressed  by  factors  0,  7,  13  and  28.  (a)  shows  rela¬ 
tive  power,  (b)  shows  signal-to-noisc  performance. 


Pulse  Compression 


In  the  compromise  between  resolution  (in  range)  and  sensitivity, 
reduction  of  the  pulse  duration  by  a  factor  a,  and  the  concomitant  widening 
of  the  bandwidth  by  the  same  factor,  reduce  sensitivity  by  20  log  a  dB. 

If  pulse  compression  is  used,  however,  compression  by  factor  b  reduces  sen¬ 
sitivity  by  only  10  log  b  dB,  i.e.  by  only  ha''.f  as  many  decibels. 

The  FPS-18  normally  has  a  pulse  duration  of  1  jisec.  We  have  installed 
a  new  pulse  transformer,  giving  a  pulse  duration  of  4  /isec.  Used  with  narrower 
bandwidth  and  without  pulse  compression,  the  change  -educes  resolution  by 
factor  4  and  increases  sensitivity  by  20  log  4  =  12  dB.  We  have  also  installed 
phase-coded  pulse  compression,  by  factor  13  or  28.  That  improves  the  resolu¬ 
tion  by  factor  13  (overall  13/4  ~  3)  and  reduces  sensitivity  by  11.4  dB 
(overall  12-11. 4  si  0).  If  that  improvement  in  resolution  had  been  achieved  by 
simply  shortening  the  pulse,  the  cost  in  sensitivity  would  have  been 
20  log  3  ~  10  dB.  Achieved  by  pulse  compression,  the  cost  in  sensitivity  is 
zero.  However,  the  improvement  in  resolution  has  limitations,  in  the  form 
of  "range  sidelobes"  down  about  15  dB  from  peak.  For  the  "extensive"  targets 
of  weather  radar,  these  sidelobes  form  a  sort  of  shelf,  down  about  15  dB  and 
extending  out  by  the  width  of  the  uncompressed  pulse:  4  ^isec  or  0.6  km. 

There  is  a  real  advantage  to  us  in  having  a  4  usee  pulse:  the  sen¬ 
sitivity  increase  of  12  dB  was  needed,  especially  for  snow.  Because  we  com¬ 
bine  data  in  range  to  reduce  fluctuations,  the  pulse  compression  affords  an 
advantage  too,  but  the  side  effects  of  range  sidelobes  are  hard  to  evaluate. 

REFERENCE:  Austin,  G.L.  and  R.W.  Fetter,  1972:  Pulse  compression 
for  meteorological  radars.  Sci.  Rept.  KW-70,  Stormy  Weather  Group,  McGill 
Univ.,  Montreal,  24  pp. 
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Standard  deviation  (cr — )  of  probability  distribution 
of  target  intensity  level,  when  k  independent  data 
are  combined  by  (a)  averaging  intensities,  (b)  averag¬ 
ing  intensity  levels,  (c)  peak  reading. 
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Dot  recording  and  peak-reading  bins 

Since  1969  our  film  records  have  consisted  of  arrays  of  discrete  dots, 
each  dot  indicating  by  its  magnitude  (both  size  and  density  change)  the  tar- 
get  intensity  -within  a  given  bin.  These  dots  are  formed  on  the  phosphor  of 
a  CRT,  and  photographed.  Discrete  values  of  dot  magnitude  are  used.  Because 
of  tube-face  flare,  the  number  of  discrete  values  must  be  kept  small:  4  or 
5  or  6.  But  it  is  because  of  tube-face  flare  that  continuous  variations  are 
hazardous,  so  that  we  '.ave  been  led  to  use  discrete  dots  and  discrete  magni¬ 
tudes. 

In  weather-radar  technique,  data  are  combined  in  groups  of  k  indepen¬ 
dent  data,  to  reduce  uncertainty  of  target  strength  attributable  to  phase 
randomness.  The  standard  deviation  of  this  uncertainty  is  given,  approximately 
at  any  rate,  by 

cr-  =  (5.57  dB)/k^ 

The  procedure  of  data-combining  affording  the  greatest  reduction  of  uncer¬ 
tainty  is  the  averaging  of  intensities.  Another,  almost  as  good,  is  to 
average  intensity  levels.  Yet  another,  which  we  use  because  of  its  simpli¬ 
city,  is  to  take  the  maximum  value  among  the  k  data.  We  use  bins  that  are 
extensive  only  in  range,  and  test  the  signal  against  a  set  of  thresholds 
spaced  uniformly  in  level,  and  report  the  highest  threshold  crossed. 

REFERENCE:  Marshall,  J.S.,  1971:  Peak  reading  and  thresholding  in 
processing  radar  weather  data.  J.  Appl.  Meteor.,  10,  1213-1233 • 
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HPPI:  Height  Plus  Plan  Indicator  (Hard  Copy) 


When  a  FASE  antenna  program  is  employed,  a  PPI  display  on  a 
persistent-phosphor  tube  face  is  not  much  use,  for  one  must  be  prepared  to 
wait  for  the  low-elevation  part  of  the  FASE  cycle.  An  improvement  can  be 
effected  by  taking  a  Polaroid  picture  of  the  low-elevation  PPI.  We  have 
improved  on  such  Polaroids  in  four  ways: 

1.  Thermofax  with  silver-salt  paper  is  better  than  Polaroid  in  that  it  can 
be  automated,  and  provides  a  bigger  picture  more  cheaply. 

2.  A  CAPPI  map  (to  range  120  nmL,  for  heights  between  5  kft  and  10  kft) 
reduces  ground  echoes  as  compared  with  a  low-elevation  PPI. 

3.  A  height  ring  around  the  map  provides  heights  in  5 -kft  intervals,  for 
the  range  interval  20  to  80  nmi. 

4.  The  information  is  fed  to  the  HPPI  display  at  a  slow  rate  over  a  period 
of  5  minutes,  so  that  the  hard  copy  can  be  generated  at  the  far  end  of  a 
telephone  line.  For  a  comprehensive  remote  display,  this  is  an  improvement 
over  our  old  system,  or  any  system  that  requires  a  few  minutes  to  build  up 
the  comprehensive  display,  followed  by  another  few  minutes  for  the  narrow- 
band  transmission.  This  narrow-band  use  of  HPPI  hard  copy  has  proved  effect 
in  a  two-year  operational  test. 

REFERENCE:  Marshall,  J.S.  and  E.H.  Ballantyne,  1970:  McGill  Weather 
Radar:  New  Research  and  Remote  Displays.  Proceedings  of  14th  Radar 
Meteorology  Conference,  Tucson,  Arizona,  p407. 


(i)  PPI,  bin  length  constant,  bin  width  proportional  to  range. 

(ii)  B-scope,  bins  of  PPI  reproduced  as  squares  of  constant  size. 

(iii)  Modified  PPI,  bin  length  equal  to  bin  width,  proportional  to  range. 

(iv)  AZEDR,  bins  of  modified  PPI  reproduced  as  squares  of  constant  size,  so  that 
ordinate  scale  is  linear  in  log  range.  (AZLOR  stands  for  AZimuth  LOg  Range.) 
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AZLOR  (Azimuth/Log  Range) 


AZLOR  in  principle:  The  radial  lines  of  a  PPI  provide  an  undistorted 
map.  For  the  diagram  in  the  lower  left,  we  have  assumed  beamwidth  1°,  one 
line  of  the  PPI  for  each  beamwidth,  and  have  drawn  every  15th  line.  We 
have  further  assumed  the  total  range  to  be  divided  into  150  bins  of  uniform 
length,  and  have  drawn  circles  at  uniform  intervals  of  15  bins.  Thus  each 
enclosure  of  the  figure  contains  225  bins,  each  having  shape  the  same  as 
the  enclosure,  but  area  less  by  factor  225.  If  each  bin  were  represented 
on  the  display  by  a  dot  (as  in  1970)  then  at  long. range  the  dots  would  be 
too  well  spaced  to  merge  into  a  half-tone  pattern,  while  at  short  range  they 
would  overlap. 

At  the  top  left  is  a  "B-scope"  display,  of  range  against  azimuth.  The 
keystone-shaped  enclosures  of  the  PPI  (a,b,c,A,B,C)  all  become  squares  on 
the  B~3cope,  whereas  on  the  undistorted  map  of  the  PPI  keystones  a,  b,  c 
were  narrower  than  the  keystone  height,  and  A,  B,  C  were  wider. 

Bin  width  is  proportional  to  range.  In  the  diagram  in  the  lower 
right,  bin  length  has  beetf  made  equal  to  width,  so  that  all  bins  are  the 
same  shape,  approximately  square.  Thus  resolution  in  range  is  everywhere 
equal  to  resolution  cross-range,  which  is  a  good  arrangement.  Now,  however, 
dots  representing  bins  would  be  very  crowded  toward  the  centre,  their  con¬ 
centration  varying  inversely  as  range  squared. 

AZLOR  equates  bin  length  to  bin  width  (as  lower-right)  but  displays 
the  dots  in  a  uniform  rectangular  array  (upper  right)  having  coordinates  of 
azimuth  and  loge  (range),  so  that  every  bin  ard  every  enclosure  of  the 
lower-right  diagram  is  reproduced  on  the  display  as  a  true  square.  The 
uniform  spacing  of  dots  means  that  the  display  cam  be  read  as  a  half-tone. 
While  gross  pattern  elements  are  seriously  distorted,  the  mapping  is  con¬ 
formal,  so  that  shape  is  preserved  for  small  pattern  elements. 

AZLJR  complements  PPI  or  CAPPI  or  HPPI:  a  person  refers  to  the 
map-type  display  for  gross  features,  transfers  to  AZLOR  for  detail.  Useful 
conformity  to  the  true  pattern  is  found  for  regions  as  large  as  r  x  r 
at  average  range  r.  22 


Display  of  Target  Intensity 


The  AZ10R  display  has  been  built  in  duplicate,  in  lower  and  upper 
sections.  In  each  section,  the  dots  are  displayed  with  four  discrete 
magnitudes  (brightness  or  density,  and  size).  The  dot  magnitudes  locate 
the  target  intensity  level  within  intervals  between  thresholds:  intervals 
of  20  aB  in  the  upper  section,  intervals  of  5  dB  within  the  20-dB  intervals 
in  the  lower  section.  They  thus  function  in  a  scale  of  four,  the  upper 
section  providing  the  fours  digit  and  the  lower  section  the  units  digit, 
to  identify  15  intervals  of  target  intensity  level,  each  of  extent  5  dBz. 


Fours  Units 

10  log  Z 

z 

R 

mm6/® 

mm/h 

♦ 

33 - 

% 

77 

5.1 

2402 

*  32 - 

x  107 

• 

72 

1.6 

1170 

31  - 

• 

67 

5.1 

570 

- 30 - 

x  106 

• 

62 

1.6 

277 

23  - 

• 

57 

5.1 

135 

•  22  - 

x  10s 

• 

52 

1.6 

66 

21 - 

% 

47 

5.1 

32 

- 20  - 

x  104 

% 

42 

1.6 

15.6 

13 - 

* 

37 

5.1 

7.6 

•  12  - 

x  103 

* 

32 

1.6 

3.7 

11  - 

* 

27 

5.1 

1.80 

10 

x  102 

* 

22 

1.6 

0.87 

03 - 

17 

5.1 

0.43 

.  02 - 

x  10 

** 

12 

1.6 

0.21 

01 

* 

7 

5.1 

0.10 

- 00 - 

-  BACKGROUND  - 


Facing:  Selected  matching  sections  of  earlier  AZLOR,  fours  above  and  units  below. 


20 


Constant-altitude  AZLOR 


Some  of  the  virtues  of  a  GP  (geometric  progression)  antenna  program 
become  evident  if  the  program  is  drawn  against  logarithmic  scales  of  height 
and  range.  For  constant-altitude  AZLOR  maps,  slices  equal  in  A  log  r  are 
selected  from  the  sequence  of  constant-elevation  maps.  The  spacing  of  the 
AZLOR  dots  ( &2r)  ^as  been  na<^e  ecIual  to  A  s,  the  spacing  between  successive 
range  lines,  so  that 

^2r  £s  2rr  1 

.  —  “  r  “S§0  -  953 

A2t  =  r/95.5 

A2  log  r  =  0.4343/95.5  -  0.00452 

Then  the  GP  factor  in  the  antenna  program  (approximately  7/6)  has  been  so 
chosen  that  A^  log  r  is  an  integral  multiple  of  &2  log  r: 

A  ^  1°6  *  **  15  A2  log  r  «*  0.678 

A  simple  reorganization  of  the  "A^"  elements  would  peim_t  an  output  of  about 
20  constant-altitude  AZLOR  maps  in  place  of  the  constant-elevation  maps  that 
come  naturally.  A  related  reorganization  of  these  "Ay  elements  would  array 
them  in  constant-range  (within  the  factor  7/6)  columns. 

Some  years  ago  we  aspired  to  achieve  this  data  processing  on  mag¬ 
netic  discs,  and  to  record  the  reorganized  data  on  magnetic  tape,  one  height 
for  each  track  on  the  tape.  We  foresaw  the  tape  record  as  highly  amenable 
to  further  specific  processing,  using  the  same  magnetic  discs.  These 
aspirations  have  not  been  achieved. 


HARPI  (Height  Azimuth  Range  Position  Indicator) 


V'e  went  to  a  4-^isec  pulse  with  two  aims:  to  use  the  pulse  uncom¬ 
pressed  for  maximum  sensitivity,  and  compressed  by  factor  i3  for  maximum 
resolution.  The  compression  factor  is  not  very  different  from  15,  the 
number  of  dots  in  a  GP  element.  Thus  we  could  obtain  about  the  same  k  by 
using  the  GP  element  as  bin  length  and  the  uncompressed  pulse  and  by  using 
the  dot-spacing  as  bin  length  and  the  compressed  pulse.  The  processed 
data  for  the  full  five-minute  cycle  can  be  displayed  on  one  HARPI  frame. 

In  this  version  of  HARPI,  the  scales  of  height  and  range  become 
logarithmic.  The  lowest  of  the  18  strips  is  for  range  15  to  ^  x  15  =  17*5  km. 
The  next  goes  from  17-5  to  ^  x  17.5  =  20.5  km,  and  so  forth,  so  that  the 
topmost  strip  is  for  range  interval  183  to  250  km.  For  the  lowest  strip, 
the  heights  range  from  1.0  to  9-0  km,  while  for  the  topmost  strip  the 
height  range  is  from  4.0  to  43-3  km.  The  hope  is  to  have  the  same  height 
scale  at  all  ranges,  say  2.0  to  18  km. 
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HEIGHT  (km)  AT  RANGE  45km 


Cloud  Photogrammetry 


For  just  a  few  hours  each  summer,  it  is  possible  to  obtain  good  stereo 
cloud  pictures.  The  following  is  the  abstract  of  a  manuscript  being  (resub¬ 
mitted  for  journal  consideration. 

Showers  observed  by  stereo  cameras  and  radar 

R.W.  Shaw  and  J.S.  Marshall 
Department  of  Meteorology,  McGill  University 

Panoramic  (rctating-lens)  cameras  with  fields  140°x55°  have  been  used 
for  stereoscopic  photography  of  clouds,  in  conjunction  with  10-cm  radar 
observations.  The  photographs  are  linear  in  the  parameters  used  in  photo¬ 
grammetry:  azimuth  and  the  tangent  of  the  angle  of  elevation. 

On  26  July  1969,  three  isolated  showers  moved  in  the  direction  (the 
same  at  all  heights)  of  the  wind.  Large  and  largely  irregular  changes  were 
noted  in  the  four-minute  interval  between  observations,  but  strong  gradients 
in  the  precipitation  observed  by  radar,  and  thin  cloud  envelopes  (from 
precipitation  to  cloud  boundary),  were  frequently  observed,  and  tended  to 
lie  in  the  general  direction  of  the  wind  and  the  late  afternoon  sun.  All 
three  shower  clouds  had  heights  of  10  km;  in  two  the  precipitation  extended 
almost  to  the  top  of  the  cloud,  but  in  the  third,  for  which  directions  of 
sun  and  wind  were  identical,  the  precipitation  top  was  at  from  5  to  7  km. 

Motion  could  be  interpreted  in  terms  of  a  mixture  ranging  from 
almost  wholly  sub-cloud  air  on  the  side  generally  toward  sun  and  wind  to 
almost  wholly  upper-level  air  for  the  anvil  tip  and  detached  anvils.  This 
indicates  that  sub-cloud  air  from  the  sunny  side  of  the  shower  was  brought 
into  the  cloud  via  an  updraught  near  the  upwind  end  of  the  shower,  and  that 
it  then  mixed  with  ambient  air  and  moved  downwind  relative  to  the 
undiluted  core  of  the  updraught. 


Facing:  26  July  1969 

Light  outline:  Visual  and  photographic  cloud. 

Heavy  outline:  Radar  isopleths,  in  steps  of  10  dBZ. 
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Radio  Lightning  Location 


▼ 


Lightning  flashes  may  be  located  by  observing  sferics  using  three 
spaced,  loop  aerial,  direction  finders.  Positions  of  events  found  in  this 
way  have  been  plotted  on  radar  maps  of  precipitation.  The  conclusions 
from  the  results  which  are  neither  as  accurate  nor  as  numerous  as  we  would 
like  are  that  the  discharges  tend  to  be  located  close  to  but  not  inside 
convective  cells  containing  the  intense  precipitation. 

The  experiment  is  being  continued  with  considerably  improved 
equipment. 

REFERENCE:  Austin,  G.L.  and  E.J.  Stansbury,  1971:  The  location 
of  lightning  and  its  relation  to  precipitation  detected  by  radar. 

J.  Atmosph.  Terr.  Phys.  22.,  841-844. 
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Example  of  slant  PPI,  including  ADA  display 
Elevation  is  8.5°  ar.d  total  range  is  55  mi. 


Microwave  Attenuation 


Eased  on  empirical  relations  among  rainfall  rate,  reflectivity,  and 
attenuation,  ADA  is  an  instrument  used  in  connection  with  our  radar  which 
provides  a  display  of  the  attenuation  due  to  rain  at  10  GHz  over  the  radar 
line  of  sight.  The  attenuation  is  indicated  as  a  deflection  trace  around 
the  perimeter  of  an  otherwise  conventional  PPI.  In  the  example  shown  on 
the  facing  page,  the  strong  echoes .toward  the  east  are  seen  to  cause  attenua¬ 
tion  in  excess  of  10  dB  (at  10  GHz)  over  an  azimuth  extent  of  about  5°. 

The  deflection  is  produced  by  an  analog  circuit  which  evaluates  an 
integral  of  reflectivity  over  range,  of  the  form 

A  =  J K  Z*  dr 

where  K  and  o<  are  empirical  constants  chosen  for  10  GHz,  and  the  limits 
of  integration  are  5  mi  and  55  mi.  The  inner  cutoff  at  5  mi  eliminates 
spurious  effects  due  to  ground  clutter. 

From  a  long  sequence  of  such  pictures,  statistics  are  compiled  on 
the  frequency  of  occurrence  of  attenuation  as  a  function  of  elevation  and 
bearing,  and  the  distribution  of  azimuth  extents  of  attenuating  regions. 

Results  from  three  summers  of  operation  show  that  the  frequency  of 
occurrence  diminishes  rapidly  as  the  elevation  increases  from  horizontal  to 
10°,  and  is  essentially  constant  from  there  on  to  20°.  5  dB  is  exceeded 
10  hr/yr  at  3°  elevation,  4  hr/yr  at  10°,  and  3-5  hr/yr  at  20°.  At  all 
elevations  10  dB  is  exceeded  about  one-fifth  as  frequently  as  5  dB. 

The  azimuth  extents  of  attenuating  regions  vary  from  1°,  the  nar¬ 
rowest  extent  measurable,  to  values  exceeding  60°.  At  low  elevations,  the 
distributions  of  azimuth  extent  are  strongly  peaked  at  the  narrow  end.  With 
increasing  elevation  the  distributions  flatten.  The  average  azimuth  extent 
of  a  5-dB  region  increases  approximately  linearly  from  10°  at  an  elevation  of 
3°  to  30°  at  an  elevation  of  20°.  For  all  elevations  the  average  10-dB 
extents  are  about  half  the  5-dB  extents. 

The  average  duration  of  a  5-dB  event  decreases  from  about  30  min  at 
near-horizontal  paths  to  less  than  10  min  at  an  elevation  of  20°.  For  all 
elevations  the  average  duration  of  a  10-dB  event  is  about  two-thirds  that  of 
the  5-dB  event. 

Papers  scheduled  for  June  publication  in  Radio  Science: 

ADA:  An  instrument  for  measuring  attenuation  due  to  rain  over  slant  paths, 
by  I. I.  Zawadzki  and  R.R.  Rogers. 

Radar-derived  statistics  on  slant-path  attenuation  at  10  GHz,  by  R.R.  Rogers. 
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A  complete  set-  of  horizontal  slices  through 
a  shower  at  the  heights  and  times  shown. 


Interpolation  and  extrapolation  techniques 
for  radar  precipitation  patterns 

The  FPS-18  radar  as  used  at  McGill  takes  5  minutes  to  record  the 
three  dimensional  structure  of  precipitation  patterns.  Specifically  for 
use  in  comparison  with  microwave-link  attenuation,  it  was  desired  to 
have  precipitation  patterns  close  to  the  surface  with  a  time  resolution 
of  about  one  minute.  Shorter  intervals  are  not  justified,  when  radar 
resolution  is  taken  into  account.  Procedures  have  been  established 
which  enable  this  coarse  sampling  time  to  be  improved  by  ir.terp  'lation 
of  patterns  from  an  earlier  time  but  higher  altitude.  Ariays  of  data 
like  the  one  shown  opposite  allow  pattern  matching  techniques  to  give 
an  'effective  pattern  fallspeed’  and  the  precipitation  pattern  to  be 
extrapolated  to  the  ground  with  a  time  resolution  improved  by  a  factor 
of  about  five.  The  sawtooth  line  (below)  gives  the  expanded  time  scale 
for  the  13  cycles  of  the  figure. 
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Achievements  in  1971-2 


1.  First  high-resolution  AZLOR  display  developed  and  put  to  use.  Film 
record  indicates  intensity  by  two  digits  of  scale  of  four  on  conformal 
maps. 

2.  First  use  of  pulse  compression  in  weather  radar  provides  high  resolution 
without  undue  sensitivity  loss. 

3.  Long  pulse  without  compression  provides  maximum  sensitivity,  recorded 
on  new  HARPI  display. 

4.  Most  severe  storms  from  summer  AZLOR  records  analyzed,  primarily  to  pro¬ 
vide  microwave-attenuation  statistics. 

5.  Storm  for  which  cloud  photographs  were  available  analyzed  to  relate 
precipitation  seen  by  radar  to  cloud  seen  by  camera. 

6.  Precipitation  measurements  by  radar  extrapolated  to  surface  (new  technique) 
and  compared  with  raingauges. 

7.  Ten  thousand  lightning  flashes  radio-located  (as  to  direction),  and  loca¬ 
tions  compared  with  radar  maps  of  precipitation. 

8.  One-hour-forecast  possibilities  explored  on  PPHI  records. 

9.  Paper  published  on  Peak  Reading  as  signal-processing  technique. 

10.  Report  published  on  interrelationship  of  showers  and  contiguous  rain. 

11.  A  thesis  was  submitted  on  autocorrelation  analysis  of  radar  maps  of  rain. 

12.  Fortnightly  seminars  for  critical  examination  of  dynamics  and  microphysics 
of  severe  storms. 

13.  Paper  accepted  on  measurement  of  snowfall  by  (CPS-9)  radar. 

14.  Paper  submitted  on  showers  observed  by  camera  and  radar. 
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Continuing  and  Future  Projects 


1.  MARSHALL  Report  and  paper  to  be  prepared  on  five-dimensional  high- 

1972 

resolution  weather  radar  system  incorporating  the  various 
McGill  technical  developments. 


2.  DRUFUCA 
AUSTIN 
BALLANTYNE 
MARSHALL 
1972 

3.  BALLANTYNE 
STAFF 

1972 

4.  BALLANTYNE 
STAFF 

1973 

5.  BALLANTYNE 
MARSHALL 
AUSTIN 
DRUFUCA 

1974 


New  techniques  of  electronic  calibration  of  weather  radar 
to  be  developed. 


Minor  improvements  in  AZLOR  display. 


Improved  version  of  AZLOR  display  to  be  built. 


Improved  amplifier  circuitry  to  be  developed,  incorporating 
as  much  as  possible  of  logarithmic  thresholds,  sensitivity  time 
control  (including  peak-reading  correction).  It  should  be 
possible  to  use  full  sensitivity  of  system  at  all  ranges. 


6.  STANSBURY 
BALLANTYNE 
STAFF 
1972  or  3 

7.  MARSHALL 
DRUFUCA 
STUDENT 
1972-3 

8.  DRUFUCA 
STUDENT 
1972  and  3 

9.  MARSHALL 
AUSTIN 
DRUFUCA 
STUDENTS 
1972  et  seq 

10.  AUSTIN 
STUDENT 
1972  et  seq 

11.  DRUFUCA 
STUDENT 
1972  et  seq 

12.  STANSBURY 
STUDENT  • 
1972  et  seq 

13-  AUSTIN 
STUDENT 
1972  et  seq 


Radio  location  of  lightning  flashes  to  be  extended  again  to 
triangulation  from  2  or  3  stations. 

5D  census  of  summer  storms  at  Montreal,  with  relevance  to 
physical  meteorology  and  application  to  microwave  attenuation. 

Attenuation  over  microwave  link:  calculations  from  radar  record 
to  be  compared  with  actual  (Bell-North  n  Research)  link. 

Thunderstorms  selected  from  summer  records  to  be  analyzed  for 
life  histories  (storm  dynamics)  and  for  storm-by-storm  approach 
to  attenuation  statistics. 

Stereo  cloud  photographs  to  be  recorded  when  weather  permits, 
and  analyzed  to  supplement  (9). 

Rainfall  and  snowfall  rates  and  amounts  to  be  determined  by 
radar,  checked  by  gauges. 

Radio  location  of  lightning  flashes  to  be  related  to  precipita¬ 
tion  pattern  by  radar. 

One-hour-forecast  techniques  to  be  developed,  introducing  storm- 
dynamic  processes.  (In  summer,  for  example,  cells  do  not 
persist  but  showers  do.)  0  . 
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Parers  In  Proceedings  of  14th  Radar  Meteorology  Conference 


1.  Comparison  of  Radar  with  Network  Gauges 

Desautels,  G.  and  X.L.S.  Gunn,  p  239 

2.  Problems  of  Snowfall  Measurement  by  Radar 

Carlson,  Paul  E.,  p  245 

3«  Motion  of  Radar  Echoes  within  Cloud  Envelope 
Shaw,  R.W.,  p  297 

4.  Finding  the  position  of  lightning  flashes 

Austin,  G.L. ,  p  313 

5.  Showers  and  Continuous  Precipitation 

Zwack,  P.  and  C.  Anderson,  p  335 

6.  McGill  Weather  Radar:  New  Research  and  Remote  Displays 

Marshall,  J.S.  and  E.H.  Ballantyne,  p  407 

7.  Radar  Weather  Performance  Enhanced  by  Pulse  Compression 

*  Fetter,  R.W.,  p  413 

6*  Autocorrelation  of  Weather  Patterns  by.  an  Incoherent  Optical  Method 
Zawadzki,  1. 1.,  p  421  $ 

9.  Decorrelation  Time  of  Weather  Radar  Signals 
Reid,  J«D>|  p  419 
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Papers  directly  dependent  on  FPS-18 


Austin,  G.L.  and  R.W.  Fetter,  1972:  Pul3e  compression  for  meteorological 
radars.  Stormy  Weather  Group  Scientific  Report  MW-70. 

Austin,  G.L.  and  E.J.  Stansbury,  197].:  The  location  of  lightning  flashes. 

J.  Atmosph.  Terr.  Phys.,  Vol.  33,  p841. 

English,  Marianne,  Charles  Warner  and  W.  Hitschfeld,  1972:  Observations 
and  theory  of  a  hail  storm.  Journal  de  Recherches  Atmospheriques 
for  the  Henri  Dessens  Memorial  issue  (submitted). 

Marshall,  J.S.,  1971:  Peak  reading  and  thresholding  in  processing  radar 
weather  data.  J.  Appl.  Meteor.,  Vol.  10,  17  pp. 

Marshall,  J.S.  and  R.W.  Shaw,  1972:  Showers  observed  by  stereo  cameras 
and  radar.  (Submitted  to  Monthly  Weather  Review.) 

Rogers,  R.R. ,  1972:  Radar-derived  statistics  on  slant-path  attenuation 
at  10  GHz.  Scheduled  for  June  publication  in  Radio  Science. 

Zawadzki,  I. I.  and  E.H.  Ballantyne,  1970:  HARPI  -  a  new  weather  radar 
display.  Q.J.  Roy.  Met.  Soc.,  £6,  144-5. 

Zawadzki,  I. I.  and  R.R.  Rogers,  1972:  ADA:  An  instrument  for  measuring 
attenuation  due  to  rain  over  slant  paths.  Scheduled  for  June 
publication  in  Radio  Science. 

Zwack,  Peter  P.,  1971:  Showers  and  continuous  precipitation:  25  July  1969. 
Stormy  Weather  Group  Scientific  Report  MW-65. 

Also  papers  from  the  14th  Radar  We ather  Conference,  1970. 
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